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Self-Consistent RPA is extended in a way so that it is compatable with a variational ansatz for 
the ground state wave function as a fermionic many-body vacuum. Employing the usual equation 
of motion technique, we arrive at extended RPA equations of the Self Consistent RPA structure. In 
principle the PauU principle is, therefore, fully respected. However, the correlation functions entering 
the RPA matrix can only be obtained from a systematic expansion in powers of some combinations 
of RPA amplitudes. We demonstrate for a model case that this expansion may converge rapidly. 
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I. INTRODUCTION 

The many body problem still being to a great deal 
an unsolved problem, new vistas are welcome. We here 
present an extension of Self- Consistent RPA-'^ which is 
based on a fermionic many body wave function being the 
vacuum to an extended RPA destruction operator. Let 
us expose the main ideas before entering the details of 
the theory and making an application to a model case. 
In this work we will consider interacting fermions, how- 
ever, the extension to the bosonic case is possible. 
Many-body theories have a clear cut hirarchy: one body, 
two body, etc. One body theories are mean field the- 
ories. One can write down a variational wave function 
for the ground state 
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where ZpJ^'' are variational parameters and \4>) — 

Y\ii<^h\vac) is an arbitrary Slater determinant built from 
an antisymmetric product of N fermion creators in the 
lowest levels, the so-called hole (h) levels below the Fermi 
energy. One can explicitly evaluate the expectation value 
of, e.g. a two body Hamiltonian 
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This 



and minimise with respect to the parameters z^^ 
leads to the well known Hartree-Fock (HF) equations. 
That this program can be pulled through hinges to a 
large extent on the fact that \Z^^'^) is the vacuum of 
quasiparticle operators ap oc exp{— Z^^'>)cpexp{Z^^^) = 

Cp + J2h^ph^h and similiar for aj^, with a|jZ(i)) = 

apl^*-^-') — 0. The scheme clearly is of the Raleigh-Ritz 
(RR) form stating that the so obtained energy is an up- 
per bound to the true ground state energy. 
Unfortunately, on the next level, i.e. the two body corre- 
lations, there does not exist such a nice closed formalism 



with RR character. As a consequence several approaches 
to two body and higher order correlation functions have 
been elaborated. Some try to be RR variational, some 
not. For example, the so-called hypernetted chain or Jas- 
trow approach puts a local, usually translational invari- 
ant, two body operator in the exponent, with u(r — r') 
the corresponding variational parameters. This leads to 
a chain of equations, the solution of which may even- 
tually converge to the exact RR variational solution of 
the corresponding energy, depending on at which level 
the chain is broken off^. The Gutzwiller ansatz'^ is tay- 
lored to the electron gas on a lattice and puts another 
two body operator in the exponent penalising double 
occupancies. However, the variation usually cannot be 
pulled through exactly. Most of the time one is satis- 
fied with the approximate solution given by the so-called 
Gutzwiller equations^. Coupled Cluster Theory (CCT) 
is an approach^ starting also with a two body opera- 
tor in the exponent (plus eventually higher order oper- 
ators) but does not lead to a RR scheme, at least not 
in its original form. Other non RR schemes are based 
on the BBGKY hirarchy of time dependent density ma- 
trices (TDDM). On this subject, there have been some 
interesting developments, recently^. In the past, several 
authors, including ourselves, tried to improve RPA as 
much as possible to overcome certain deficiencies for the 
calculation of two body correlation functions. It is, e.g., 
well known that standard RPA (s-RPA) yields two body 
correlations functions which violate the Pauli principle 
due to the so-called quasi-boson approximation inherent 
to s-RPA (we define s-RPA as being the small amplitude 
limit of TDHF implying that the exchange term is in- 
cluded). Usually this leads to a quite strong overbinding 
of the ground state energy whereas excitation energies 
may be more correct. We here will formulate the RPA 
as a diagonalisation problem as is most appropriate for 
finite Fermi systems, since a good part of the spectrum 
consists of discrete levels. In this way one usually writes 
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the RPA excitation operator aai 



(3) 



with k = {p, h) the (particle, hole) levels. The excited 
state is then written as 



(4) 



where |0) is the ground state so that it is the vacuum to 
the destructors Qi^, i.e. 



(5) 



Since this scheme cannot be pulled through in general, in 
s-RPA the so-called quasi boson approximation is intro- 
duced replacing the fermion pair operators by ideal Bose 
operators^. 
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In this way, ([3]) boils down to a Bogoliubov transforma- 
tion among ideal bosons with the usual HFB vacuum^ 
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with a\ap\HF) Bph\HF) ~ 0. Since a fermion 
Hamiltonian can be expanded via Holstein-Primakoff or 
Belyaev-Zelevinsky boson expansion in powers of the bo- 
son operators^, the RPA ends up with a standard HFB 
approach for interacting bosons. The scheme approaches, 
of course, a RR variational procedure only in as much as 
the boson expansion of the Hamiltonian converges. 
In the past, we used the RPA operator ^ with the 
fermion pair operators, supposing that the killing con- 
dition Qy\Q) — existed. However, a corresponding 
vacuum can only be constructed in some very simplified 
model cases (see, e.g,— ) and, thus, in general, the hypoth- 
esis of the existence of a vacuum remained an approxi- 
mation what perturbed the internal consistency of the 
so-called Self Consistent RPA (SCRPA) approacbii^i^. 
In this work we will adopt the strategy to keep in ^ 
instead of the boson operators the original fermion pair 
operators. We then will see that this ground state still 
is vacuum to a generalised RPA operator containing the 
standard operator ([3]) plus, in addition, specific two body 
operators correcting for the Pauli principle. With that 
RPA operator we will set up the familiar equation of 
motion (EOM) method^. The problem will be to evalu- 
ate the expectation values of the commutators, since the 
new RPA operator is non linear, i.e. not of the form of 
a canonical transformation. Therefore, an effective sys- 
tematic expansion scheme will be set up for these expec- 
tation values. The procedure will be applied to the very 
schematic two level Lipkin model where it is shown that 
the expansion converges fast, even for very small particle 
numbers. 



II. GENERAL THEORY 

As mentioned in the introduction, we will consider the 
following variational ansatz for the many body ground 
state 
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where Zp^^^^^^y^^ are the variational parameters. In gen- 
eral they can be complex. In principle we can add a one 
body operator as in ([l} to the exponent and vary z*-^-* 
and z*-^-* simultaneously. In order to keep the formalism 
as transparent as possible, we will refrain from this com- 
plication here and suppose that we work in the standard 
HF basis. The first task will be to find an RPA operator 
(5+ which applied to \Z^'^^) gives the excited states of the 
system: \v) = Q+jZ*^^^) and whose hermitian conjugate 
kills the ground state acting, therefore, as a vacuum 



Q.|^('^) = 



V V. 
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We make the following ansatz for the destruction 
operatopiS 
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We indeed verify that kills the ground state 
The amplitudes A, Y can be normalised according to 
(0|[Qi,,(5+]|0) = 1. Evidently contains the usual 
RPA operator ^ plus extra specific 2-body operators. 
These extra terms serve to fullfil the Pauli principle. 
They make, however, the calculus rather heavy. In the 
following, we shortly outline our approach and then give 
an application to a schematic model. 
In order to determine the amplitudes A, Y, we proceed 
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as usual and minimisethe following mean excitation en- 
ergy, equivalent to the minimisation of an (symmetrised) 
energy weighted sum rule^, 



{0\[[Q,,H],Qt] 
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Using in this expression the standard RPA operator 
= -^ph'^p'^h — Ypf^a\ap, we obtain from the minimi- 
sation of (|16l) with respect to the amplitudes X and Y the 
usual structure of the (Self Consistent) RPA equations. 
However, with ([5]) the RPA equations become more com- 
plicated. Among other things, this stems from the fact 
that the norm matrix is now non diagonal, i.e. the struc- 
ture of the equations is as follows 



^11 Ri2 \ f ^ 
i?2i R22 \Y 



Nil N12 
N21 N22 
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with obvious definitions of the matrices R and N . They 
are hermitian matrices by construction, i.e. R = i?"*" and 
N = N^. In order to arrive at the usual structure, we 
write N = N^^^NqN^^"^ with 



No 



1 

-1 



and obtain in short hand notation 



RX - n.Nox 



(18) 
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with R = N-'^/^RN-^/^ and x = N'^''^X- This equa- 
tion now has globally the same mathematical structure 
as standard RPA or SCRPA (details are different, of 
course). In particular, we have for the normalisation 
X^ — Y^ = 1. Of course, the task is not finished. The dif- 
ficulty resides in how to express the correlation functions 
in R and N by the amplitudes X,Y , since the form ([12]) 
of the RPA operator, as already mentioned, does not con- 
stitute a canonical transformation among operators and 
can, therefore, not be inverted in a straightforward man- 
ner. We will set up an expansion scheme which can be 
paraphrased as follows: Suppose we did not know how to 
invert the Bogoliubov transformation from ordinary to 
quasiparticles in standard BCS theory. How to calculate 
expectation values in the BCS state nontheless? Among 
other ways, this can go as follows. As is well known, the 
BCS state can be written as a coherent state: \BCS) = 
sxp\^ Zk(^^.(^_fJ\\vac) . We want to calculate, e.g. the pair- 
ing tensor Kk = {BCS\c-kCk\BCS) / {BCS\BCS) . We 
have Ck\BCS) — ZkC^_^\BCS). Multiplying this relation 
with c-k and repeating the operation leads to the well 
known BCS expression Kk — UkVk where we supposed for 
simplicity Zk to be real, i.e. {c-kCk) = (c|.c^fe), and used 
Zk — Vk/uk with u1+vl — 1. Using an analogous method 
to calculate the correlation functions in the RPA matrix 
leads to a series which, unfortunately, does not break off 
after a low number of terms unless one deals with very 
few particle states. As the example with the BCS state 



shows however, the series contains nontrivial terms and 
may, thus, converge rapidly. If convergence is achieved, 
Pauli principle is respected fully, since the theory is based 
on a wave function. Of course, the final success will de- 
pend on that supposed rate of convergence. This ends 
the rough outline of the theory. Let us now make an ap- 
plication to the Lipkin model which often has served in 
the past as a testing ground for many body theories. We 
will see how in that model the series expansion works and 
that in this case the convergence of the series is, indeed, 
fast and performant. 



III. APLICATION TO A SCHEMATIC MODEL 

In this section we want to apply our theory to the 
schematic Lipkin modelii which is an exactly solvable 
two level model with separable monopole-monopole in- 
teraction. The Hamiltonian is given by^ (without loss of 
generality we put the intershell spacing e = 1) 
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H = Ja~-{ji + j^ 
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with V the coupling constant and 
2 
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J 
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where fl is the degeneracy of one shell. We take ft = N 
with N the particle number, i.e. the lower shell is filled 
for V = 0. The SU2 commutation relations are 



2 Jo 



Jo, J± 



± J±. 
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The variational ground state corresponding to ^ is given 
by 



with 



\Z)^e^\HF) 



zJj^J+. 



(23) 
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We verify that this correlated state is vacuum to the fol- 
lowing destruction operator— 



Q = XJ^ - YJ^ 
with z = Y/ [NX) and 

J+ = ( 1 - ■qJa] J+, 



N 
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We note that the non-standard correction term is of order 
1/A^. In this model the matrix elements of eq. pT]) are 
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given by 



Rii — 



Rl2 

R21 



\H, J\ 



J-,H 



J4 



[J+,H],J+ 



R22 = -(([j+,[iI,J-]]) + ([[j+,i/],J-])) (27) 
and the elements of the norm matrix, 



N12 



> = -2(io> 



J-,J- 



J 



N21 



J+,J+ 



N22 = {[J+.J-]) 

= 2(io) - 277 [2(Jo) + 2(j2) - (J+J_) 
+772|2(Jo) + 4(J2) + 2(J3) - 3(J+J_) 

-2(J+Joi->l. (28) 



To calculate the double commutators in ([27|) . we use the 
Jakobi identity [A, [H, B]] - [[A, H],B]^ [H, [A, B]], and 
obtain 



Rii 
R12 
R22 



-2{Jo)+V[{Jl) + {Jl) 

R21 = -V^|2(Jo) + 4(J2) _ 2(J+J_) - 77 [3(Jo) + 8(J2) + 4(J3) - 7{J+,L) - 6{J+JoJ-) 
-2(Jo) + V ((J2 ) + (J2 )) + [4(Jo) + 4(j2) + 2(J+ J_) + V ((P) + (J^ ) 
-7?2|2( Jo) + 4(J2) + 2{J^) - 3( J+J_) - 2{J+Jo.L) 



-V 



Y + (JD) + 8 ((-^oJi) + (JlJo)) + {{JpD + (JlJ^)) 2 ((4 J_) + (J. 



(29) 



Diagonalising the norm matrix we arrive at an equation 
which has the sturcture (fT9|) . i.e. the same structure as 
the standard RPA or SCRPA equations. 
At this point we shortly would like to explain how one can 
calculate the correlators appearing in the RPA matrix 
within our formalism following in spirit the example given 
in section II for the calculation of the pairing tensor. Let 
us demonstrate this on the example of {Z\Jo\Z) / {Z\Z) 
to lowest non trivial order. We have 



Thus, 



Jo\Z) 

J-\z) 

J-J-\Z) 



2z{N -\)J+\Z) + 
2z{N - 1) 



-2Jo + J+J- 



\Z) 



2z[N-l) -2Jq + 2z{N -1)J+J+ \Z) + 



Now, we calculate the average value in the state | Z) and 
use {Z\J_L\Z) = {Z\J+J+\Z) 

N - - 

{Z\Jo\Z) = --{Z\Z)+2z{Z\J+J+\Z) 



{Z\J+J+\Z) 



-Az{N - 1) 
l-iz'^iN - 1)' 



:{Z\Jo\Z). 



{Z\Jo\Z) 



N 



l-4z2(7V-l)2 



{Z\Z) 2 l~Az^{N ~1){N -3)' ^^^^ 

With this it becomes evident how to push the expansion 
to higher order and how to do analogous expansions for 
all other correlators. A more complete set of correlation 
functions is given in the Appendix (f51) . We now will solve 
eq (fT7|) for different particle numbers. For = 2, eq (fT7|) 
reads 

2(l-4z2 + 4yz) -y(l + 4z2) 
-V{l + Az^) 2{l + 2Vz) 

2(1- 4z2) 4z \ / a: 

4z -2 [ Y 



(31) 



where we divided by the common normalisation factor 
(1 + 4z^). Let us begin with the diagonalisation of the 
norm matrix 

2(1- 4z2) 4z 
4z -2 



n = 



(32) 



Its eigenvalues are 

e+ = 2Vl + 4z4-4z2 >0, 
e_ 



-2Vl + 4z4-4z2 < 



(33) 
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and the corresponding unitary matrix is given by 
with (det ([/) = 1) 



u V 
— V u 



(1 - 2z2) + x/TTI? 



4z2 + [1 - 2z2 + VI + 4z4J ' 
2z 

4z2 + [i _ 2z2 - VI + 4z4J^ 



(34) 



(35) 



The norm matrix (|32|) can be written as n = it}/'^Nqtt}/'^ 
with 



and vice versa 




(40) 



Therefore, the RPA operator can be written as 

Qu = XJ^-YJ+ 
= XJ1+YJ2 



with 



J2 



V^ V 



uJ- — vJj, 



vJ- + uJ+ 



(41) 



We thus obtain rx ~ ^uNqx oi' exphcitly 



rii ri2 
hi f22 



1 

-1 



(36) 



with X'^-Y'^ = 1 and 



rii 



ri2 = '~2i = ^ ("^ ^ ''12 + ^"""^ (''ii ~ ^22) 
?^22 = ^ (w^r22 + v'^rii + 2uvri2) . 

Multiplying first equation of ([55]) with F and second with 
X we obtain after adding and substracting both equa- 
tions 



{hi + r22) XY + fi2 + 







{hi - h2) XY + f 12 - r2 j = 2Vt^XY. (37) 

For the solution of this equation, we need to express X 
and Y by X, F. This is given below in dSH). With 
z — Y/{NX) one then obtains 



z± 



V 



1 



2 1± vr+y2 



±\/l + V2. 



(38) 



It can readily be checked that this solution coincides with 
the exact one. Please notice that the equality fn = f22 
only holds at the solution Therefore, only after 

convergence to the solution we obtain the usual RPA 
structure with eigenvalues which come in pairs ±|rJiy|. 
This also holds for the other cases with N > 2 considered 
below. We surmise that this is a general property of our 
theory. 

The new amplitudes X , Y in terms oi X , Y are given by 



(39) 



The solution for iV = 4, iV = 6 and TV = 10 can be 
obtained with the correlation functions given in the ap- 
pendix . The results are shown in Figs [1] [2l IH and 
gl We see that SCRPA with 77 ^ (SCRPA(r;)) system- 
atically and substantially improves over the rj — case 
(SCRPA (0)) in the range x < 1. This fact shows that 
the extra correlations introduced in the excitation opera- 
tor (jl2p improving the Pauli principle play an important 
role . In the range x > li ground and excited states be- 
come strongly mixed and SCRPA in the 'spherical' basis 
deviates more and more from the exact solution. The 
good agreement of SCRPA (77) with the exact solution 
also demonstrates that the convergence of our expansion 
scheme is fast. In future work it will be very interesting 
to apply our theory also to the symmetry broken regime. 



IV. CONCLUSIONS 

In this work we have set up a Self-Consistent RPA 
(SCRPA) theory which is based on a fermionic many 
body vacuum. That is, we constructed a killing oper- 
ator of the correlated vacuum |0) ^ ea;p[Z(2)]|ffF), with 
Z(2) a two body operator, see eq ([5]), (flUl) which also 
is the starting point of Coupled Cluster theory^. This 
killing operator consists of the usual ph RPA operator 
with, additionally, specific two body terms correcting for 
the Pauli principle. With this extended RPA operator, 
we develope the standard Equation of Motion technique 
and arrive at SCRPA equations of the usual mathemat- 
ical structure. Since the extended RPA operator is not 
of the type of a canonical transformation between opera- 
tors, the evaluation of the correlation functions contained 
in the RPA matrix in terms of the RPA amplitudes is not 
evident. We proposed a systematic expansion scheme in 
terms of higher and higher correlations functions (there 
may be room for improving this scheme) and applied it 
to the two level Lipkin model as a testing ground. It 
is found that the additional terms in the RPA operator 
systematically and substantially improve the results. We 
want to stress that once convergence is achieved with 
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the mentioned expansion, the Pauh principle is fully re- 
spected, since our method is based on a fermionic many 
body ground state which is the vacuum to an (extended) 
RPA operator. In this way SCRPA which in the past 
suffered from the fact that it was based on a non existing 
vacuum is now put on grounds with more internal consis- 
tency. However, the practical success will depend on the 
convergence rate of our expansion scheme in more realis- 
tic situations. This is similar to the convergence rate of 
the hirarchy of hypernetted chain equations with Jastrow 
theory. In this work we only investigated the symmetry 
unbroken phase of the Lipkin model. The symmetry bro- 
ken phase as well as other models for strongly interacting 
fermions shall be investigated in the future. On the for- 



mal level some questions still need to be considered also. 
The most important one is whether the Goldstone theo- 
rem when working in a symmetry broken phase still holds 
as in s-RPA. In particular this concerns the fuUfillment 
of conservation laws, the appearence of a zero mode and 
the respect of Ward identities. 
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Appendix A: Correlation functions 
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In this Appendix, we list correlation functions to lowest 
non-trivial order as they are needed for the RPA calcu- 
lation 
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FIG. 1: First excited state energy with s-RPA, SCRPA(O), 
SCRPA(77), and exact solution for = 4 as a function of 
X = V{N-1). 




FIG. 2: Same as Fig. [T]but for A = 6. 
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FIG. 3: Same as Fig. [T]but for = 10. 
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FIG. 4: Same as Fig. [T]but for = 20. 



